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ABSTRACT Using DFT/TDDFT methods, the excited-state lifetimes of Ru(II) polypyridyl complexes were 
computed accurately and the reason of Ru(ID) polypyridyl complexes with long excited-state lifetimes was 
explained by the electron-transfer distances and HOMO-LUMO gaps. Finally, the photovoltaic conversion 
efficiencies of complexes were predicted using DFT and docking methods. This work has provided methods of 
predicting the excited-state lifetimes and photovoltaic conversion efficiencies of Ru(II) polypyridyl complexes. 
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1 INTRODUCTION 


Owing to the excellent photophysical and photochemical properties, Ru(II) polypyridyl complexes have 
attracted considerable attention in many important areas such as photochemical conversion of solar energy, 
DNA-photocleavage reagent, molecular ‘‘Light Switch’’, etc!!, For example, Ru(II) polypyridyl complexes 
can effectively cleave DNA under light irradiation and the DNA-damage reason is attributed to the produced 
singlet oxygen (!O2)P!, usually related to the quantum yield! ?!, Besides, many Ru(II) polypyridyl complexes 


can emit strong fluorescence and their excited-state lifetimes play an important role in quantum yield. Hence, 
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studies on excited-state lifetimes have been a hot topic and are explored widely in experiments. Research 
results'®'© show that the emission intensity relates to the excited-state lifetimes of complexes, usually 
dominated by metal-centered (MC) states and metal-to-ligand charge (MLCT) states. However, these states are 
experimentally determined with difficulty. If MC and MLCT states of complexes can be obtained theoretically, 
the results will be of great significance for predicting the emission intensity and directing the design of novel 
luminescence complexes. 

In this work, we select the experimentally reported Ru(II) polypyridyl complexes!!! 1~5, i.e., 
[Ru(tpy-PhCH3)2]**, [Ru(tpy)2]**, | [Ru(H2pbbzim)2]**, [Ru(tpy-HImzphen)]** and [(tpy-PhCH3)Ru(tpy- 
HImzphen)]?* (Structural diagrams are shown in Fig. 1) to perform a theoretical study using density functional 
theory (DFT) and time-dependent theory (TDDFT) methods'?-!, In addition, as is well-known, many Ru(II) 
polypyridyl complexes can be widely used in dye-sensitized solar cells (DSSCs). To predict the photovoltaic 
conversion properties of these complexes, the excited-state properties of these complexes were also explored" 
Ul. We hope that this work laies a theoretical foundation for designing and synthesizing novel complexes with 


better luminescence properties. 


2 THEORY AND COMPUTATIONAL METHODS 


2.1 Computations of the excited-state lifetimes 

Full geometry optimization of the studied complexes 1~~5 in the ground state was carried out using the 
restricted B3LYP method with the LanL2DZ basis set'*: !?! for Ru atom and with the 6-31G(d) basis set for the 
other atoms. For the obtained structures, frequency calculations were also performed using the same method to 
verify that the optimized structure was an energy minimum. To obtain energies of MLCT and MC states, 240 
singlet-excited-state energies of these complexes were calculated based on the optimized ground geometries 
with TDDFT method at the same level of theory. 

Many Ru(Il) complexes can be excited under light irradiation and their excited-state lifetimes are 
governed by the nonradiative decay rate constant knr, given byl* ?! 

knr = Kn kar (1) 

where the overall radiationless decay knr is the sum of kn? and kn’. Specifically, kn° leads directly from the 
MLCT state to the ground state, whereas knr' is related to a thermally activated process that takes into account a 


surface-crossing to a low-lying metal-centered (MC) level, so it depends on the energy gap AE between MLCT 
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and MC states. 
The excited-state lifetimes (r) of complexes can be obtained by the reciprocals of Eq. 1, expressed as 
follows: 
T= lu. ka) (2) 
The rate constant k (ks? and ks) can be obtained by Eq. 3. 
k = Aexp[-E(RT)] (3) 
where A is the frequency factor, T is the temperature, and E. is the free activation energy. A can be obtained by 


the following equation: 


A=2ad’L T (4) 
7ü 


where d is the diameter of the complex, L is the Avogadro constant, T is the temperature, and M is the molar 
mass of the complex. The excited-state lifetimes were obtained via Eqs. 27-4. 
2.2 Computations of photovoltaic conversion properties 

The optimized structures of these complexes were docked with two TiO» molecules using the Dock6.0 
program!*°!, The box size, grid space, energy cutoff distance, and maximum orientation were set as 20, 0.3, 
9999 A, and 200,000, respectively. The other parameters used in docking were default. There are rigid docking 
and flexible docking, which represent two quite different docking approaches. In rigid docking, which is based 
on receptor spheres generated by the SPHGEN module as well as the heavy atom centers of the ligand, the 
ligand is docked rigidly to the receptor, whereas in flexible docking, which is based on the 
anchor-and-growing algorithm, the ligand is docked flexibly to the receptor. Here, using the rigid docking 
method, the docking models (complex-TiO?) were obtained. To accurately compute the electron-transfer 
properties between Ru(II) complexes and TiO», the optimizations of the obtained docking models in ground 
state were further carried out using CAM-B3LYP method and with the LanL2DZ basis set!!*!°! for Ru atom 
and with the 6-31G(d) basis set for the other atoms. The optimized results are given in Fig. 2. Meanwhile, the 
frequency calculation was also performed in order to verify the optimized docking model to be an energy 
minimum. In addition, to obtain the excited properties of the docking models, the docking models in the lowest 
triplet/singlet states were also optimized at the same level of theory. 


All calculations were preformed with the Gaussian09 program-packageP". 


3 RESULTS AND DISCUSSION 


3.1 Evaluation of computational accuracy 

Complex 4 was selected for optimization at the level of B3LYP/LanL2DZ + 6-31G(d), since the crystal 
structure of complex 4 has been determined. The calculated results of complex 4 and corresponding X-ray data 
are listed in Table 1. Comparing the calculated geometrical parameters of complex 4 with corresponding X-ray 
data, we can clearly see that the calculated results using the DFT method are in a satisfying agreement with the 
corresponding X-ray datal!!, 

In addition, complexes 4 and 5 were selected for calculations of electronic absorption spectra, since 
electronic absorption spectra of complexes 4 and 5 have been determined experimently. The computed 
electronic absorption spectra of 4 and 5 are given in Fig. 3. We can see that the computed maximum bands of 
complexes 4 and 5 are at 498.2 nm (f= 0.617) and 502.8 nm (f = 0.759), respectively, in good accordance with 
experimental results!!! 492 and 502 nm, respectively. From Fig. 3, we can also see that the simulated 
absorption spectra of 4 and 5 agree well with experimental results whither from absorption-spectra shapes or 
positions of bands. This shows further that the computed results are believable using such a DFT method. 

3.2 Excited-state lifetimes of complexes 

The energies of MLCT and MC states of Ru(II) polypyridyl complexes play an important role in the 
excited-state lifetimes. To obtain the excited-state lifetimes, absorption spectra of complexes 1~~5 were 
computed using TDDFT method and the computed energies of MLCT and MC states are given in Table 2 and 
the computed excited-state lifetimes t via Eqs. 2~4 and corresponding experimental results are given in Table 
3. We can see that the computed t of complexes 1~5 are 4599216.6, 6368.4, 202.6, 66222084.1 and 
8752385.3 ns, deviating from the experimental results very greatly. Nevertheless, the general tends are in 
accordance with the experimental results!!!l, i.e., 7(4) > (5) > «(1) > «(2) > «(3). The main reason of error may 
be that the rate constant k in Eq. 3 was obtained by classical collision theory and light irradiation is not taken 
into account. Based on the idea, Eq. 3 is revised further and expressed as follows: 

k = AA’exp[-E(RT)] (5) 

Complex 1 was selected as an example. The experimental result t for 5.0 ns was substituted in Eq. 2 and 
A' for 920037.6 was obtained via Eqs.1~5. Using the obtained A’, the revised t values of complexes 2~~5 
were computed being 0.07, 0, 71.9 and 9.5 ns, more close to the experimental results!!! of 0.25, 0, 55.5 and 
10.2 ns. Therefore, the frequency factor A' for 920037.6 may be used for computations of excited-state Ru(ID 


polypyridyl complexes and the excited-state lifetimes of Ru(II) polypyridyl complexes can be predicted 
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accurately using the frequency factor. 
3.3 Molecular orbital analysis 

Frontier molecular orbitals play an important role in exploring relation between molecular structure and 
excited-state lifetimes. So, it is necessary to analyze and discuss enengies and components of frontier 
molecular orbitals. The highest occupied molecular orbitals (HOMO), the lowest unoccupied molecular 
orbitals (LUMO) and HOMO-LUMO gaps were computed and listed in Table 4, and the corresponding 
stereocontour plots are shown in Fig. 4. 

From Fig. 4, we can see that the "electron cloud" of the HOMO in ground states is mainly distributed on 
the end of one ligand, whereas that of the LUMO on another ligand for complexes 4 and 5. This shows that 
electrons transfer from the end of one ligand to another ligand when electrons are excited, leading to the 
electron-transfer distances to be far away. On the contrary, electrons transfer from one ligand to another ligand 
or from Ru atom to ligands when electrons are excited for complexes 1~~3, leading to the electron-transfer 
distances to be short relative to complexes 4 and 5. In addition, it can be seen from Table 4 that the computed 
HOMO-LUMO gaps of complexes 1~5 are 3.222, 3.603, 3.094, 1.087 and 1.190 eV, respectively, showing 
that complexes 4 and 5 have less HOMO-LUMO gaps relative to 1~~3 and they can be excited easily. 

From the above analysis, we can see that complexes 4 and 5 with long excited-state lifetimes have two 
characteristics: (1) less HOMO-LUMO gaps; (2) long electron-transfer distances. The two characteristics may 
be the reason of excited-state lifetimes of complexes 4 and 5 longer than those of complexes 1~3. 

3.4 Photovoltaic conversion property 

The calculated net charges of two TiO» molecules in the docking models of complexes 1~~5 in the ground 
states and in the lowest triplet states are listed in Table 5. Since the photovoltaic conversion efficiencies of 
complexes 1~5 were not determined experimentally, for comparison, complexes?? 3a and 3b with 
photovoltaic conversion efficiencies in experiments were also computed and the computed net charges of two 
TiO» molecules in the their docking models in the ground states and in the lowest triplet/singlet states were 
also listed in Table 5. For complexes 3a and 3b, the calculated net charges of two TiO» molecules in the 
docking models in the ground states are —0.3079 and —0.2098 |e|, and those in the lowest singlet states are 
—0.3356 and —0.2305 |e|, respectively. Their gaps (€) are 0.0277 and 0.0207 |e|, showing that greater charges 
transfer from complex 3a to TiO» molecules relative to complex 3b and the photovoltaic conversion efficiency 
of 3a should be better than that of 3b, inconsistent with the experimental result"?!, Similarly, the calculated 


gaps of net charges of two TiO» in docking models between in the ground states and the lowest triplet states for 
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3a and 3b are 0.0217 and 0.0499 |e], respectively, showing that less charges transfer from complex 3a to TiO» 
molecules relative to complex 3b and the photovoltaic conversion efficiency of 3b should be better than that of 
3a, in accordance with the experimental result??! Therefore, the net charges in the docking models of 
complexes 1~5 in the lowest triplet states were computed below. 

From Table 5, we can see that the gaps (s) of net charges on two TiO»in docking models between in the 
ground states and in the lowest triplet states are 0.0127, 0.0237, 0.1256, 0.0022 and 0.0287 |e|, respectively. We 
predict that the order of photovoltaic conversion efficiencies of complexes 1~5 should be3 > 5» 2» 1» 4. 
Such a result is inconsistent with the excited-state lifetimes of complexes 1~5. This further shows that 


photovoltaic conversion efficiencies of complexes are not related with their excited-state lifetimes'?l, 


4 CONCLUSION 


Theoretical studies on the excited-state lifetimes and photovoltaic conversion efficiencies of Ru(II) 
polypyridyl complexes 1~~5 have been carried out using the DFT/TDDFT and docking methods, leading to the 
following results: (1) the calculational method of excited-state lifetimes for Ru(II) polypyridyl complexes 1~~5 
was explored and the excited-state lifetimes can be predicted accurately. (2) The reason of excited-state 
lifetimes of complexes 4 and 5 longer than those of 1—3 was explained by the electron-transfer distances and 


HOMO-LUMO gaps. (3) The photovoltaic conversion efficiencies of complexes 1~~5 were predicted. 
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Fig.1. Structural diagrams of complexes 1~5 and atom labels 


Fig.2. Optimized docking models of complexes 1—5 with two TiO2 molecules 
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Fig. 3. Simulated absorption spectra of complexes 4 and 5 
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Fig. 4. Molecular orbitals of complexes 1~5 in the ground states 


Table 1. Computed Selective Bond Lengths (A) and Bond Angles (°) of 


Complex 4 and Corresponding Experimental Data"!!! 

Bond Dist. Expt. Angle (°) Expt. 
Ru(1)-N(2) 2.1303 2.049 N(2)-Ru(1)-NG) 78.04 78.84 
Ru(1)-N(3) 2.0192 1.966 N(3)-Ru(1)-N(4) 78.05 79.31 
Ru(1)-N(4) 2.1308 2.054 N(5)-Ru(1)-N(6) 7131 78.03 
Ru(1)-N(5) 2.1404 2.069 N(6)-Ru(1)-N(7) 7138 78.84 
Ru(1)-N(6) 2.0517 2.016 N(5)-Ru(1)-NQ) 92.79 92.43 
Ru(1)-N(7) 2.1397 2.066 N(7)-Ru(1)-N(4) 92.73 92.29 


Table2. Calculated Related Excitation Energies (AE/eV), Oscillator Strengths and 
Main Orbital Transition Contributions of Complexes 1~5 Using the TDDFT Method 


Major contribution AE/eV A/nm f Character 


1 H-1—L-5 (45.1%) 3.7572 329.9 0.136 d—2a* (MLCT) 


H-L+6 (48.1%) 


H-1L-*11(37.296) 4.1062 301.9 0.023 d—d* (MC) 

2 H—>L+2 (79.4%) 3.0899 401.3 0.112 d—* (MLCT) 
H-2—L (27.3%) 2.9153 425.3 0.043 d—d* (MC) 
H-1—L+1 (27.3%) d—d* (MC) 

3 H-1—L-43 (100%) 3.0121 411.6 0.013 d—2a* (MLCT) 
H-2—L-^1 (45.9%) 2.9203 424.6 0.073 d—d* (MC) 


H-1—>L (45.9%) 


4 H-7—L-^1 (100%) 2.8757 431.4 0.033 d—* (MLCT) 
H-8—L+2 (60.8%) 3.2919 376.6 0.011 d—d* (MC) 

5 H-7—L22 (90.4%) 3.0453 407.1 0.136 d—n* (MLCT) 
H-4—L (44.8%) 2.6824 462.2 d—d* (MC) 


Table 3. Calculated Excited-state Lifetimes (rt, ns), Corrected Excited-state Lifetimes (7*P, ns) 


and the Experimental Values"! of Complexes 1~5 


1 2 3 4 5 
T 4599216.6 6368.4 202.6 66222084. 1 8752385.3 
Expt. 5.0 0.25 0 39:0 
pP 5.0 0.07 0 71.9 9.5 


Table 4. Computed Energies (eV) of HOMO, LUMO and HOMO-LUMO Gaps (e) 


Enowo Etuvo E 
1 —10.308 —1.086 3.222 
2 —11.186 —1.583 3.603 
3 —10.278 —1.184 3.094 
4 —8.130 —1.043 1.087 
5 —8.120 —6.930 1.190 


Table 5. Computed Net Charges (|e|) of Two TiO: in Docking Models in the Ground States 


and in the Lowest Singlet/triplet States and Gaps (£) for Complexes 1—-5, 3a and 3b 


Comp. 3a 3b 1 2 3 4 5 
Charge? -0.3079 -0.2098 -0.2601 -0.1532 -0.3309 -0.4061 -0.1569 
Charge” -0.3296 -0.2597 0.2728 0.1769 —0.4565 —0.4083 0.1856 
Charge* -0.3356 -0.2305 

E 0.0217 0.0499 0.0127 0.0237 0.1256 0.0022 0.0287 
e 0.0277 0.0207 


a express the net charges of two TiO» in docking models in the ground states; > express the net charges of two TiO» in docking models in the lowest triplet 


states; * express the net charges of two TiO» in docking models in the lowest singlet states; 4 express the gaps of net charges on two TiO» in docking models 


between in the ground states and in the lowest triplet states; * express the gaps of net charges on two TiO» in docking models between in the ground states 


and in the lowest singlet states 
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The excited-state properties of Ru(II) polypyridyl complexes were explored and their excited-state lifetimes 
and photovoltaic conversion efficiencies were predicted. This work has provided methods of predicting the 


excited-state lifetimes and photovoltaic conversion efficiencies of Ru(II) polypyridyl complexes. 
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